16-1 


DROPLET  SIZE  DISTRIBUTION  AND  LIQUID  VOLUME  CONCENTRATION 
IN  A WATER  SPRAY:  PREDICTIONS  AND  MEASUREMENTS 

Gabriel  P.  Pica 
Inscituto  Superior  Tecnico 
Dep.  Engenharia  Mecanica 
1096  Lisbon 
Portugal 

AD-P003  141 


SUMMARY 

*" s’  A water  spray,  from  a twin-fluid  atomizer,  was  studied  at  different  air  temperature  and  at  atmospheric 
pressure. 


The  Sauter  Mean  Diameter,  D^,  and  the  droplet  volume  distribution  were  measured  at  four  different 
distances  from  the  injector. 

An  optical  method  was  used  to  measure  droplet  sizes. 

Predictions  of  the  Sauter  Mean  Diameter,  liquid  volume  concentration  and  droplet  size  distribution  were 
also  evaluated.  The  influence  of  the  air  velocity,  air  temperature  and  of  the  water  flow  rate  in  the  spray 
S.M.D.  and  in  the  liquid  volume  concentration  at  different  distances  from  the  injector  have  been  measured. 
The  predictions  showed  good  agreement  with  experimental  results.  , 
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Droplet  surface  area 
Specific  heat  of  droplet  liquid 
Liquid  volume  concentration 
Spray  Sauter  Mean  Diameter  (S.M.D.) 

Heat  transfer  coefficient 
Latent  heat  of  vaporization 
Droplet  mass 
Prandtl  number 
Droplet  radius 
Droplet  Reynolds  number,  Re 

- K 

Tm  - Air  temperature 

- Droplet  temperature 

ure£  - Relative  velocity  of  the  droplet  to  the  mean  stream 


pureltDd 


p,U  - Gas  properties  evaluated  a film  temperature 
p„  - Liquid  density 


(- 


1.  INTRODUCTION 

Knowing  the  behavior  of  an  evaporating  liquid  droplet  is  essential  to  the  understanding  and  prediction 
of  the  performance  of  the  spray. 

The  behavior  of  the  liquid  fuel  spray  plays  an  important  role  in  the  combustion  efficiency  and  in  the 
production  of  the  pollutants  in  combustion  chambers. 

Theoretical  studies  of  evaporating  sprays  have  been  carried  out  by  several  authors  | 1,2, 3,4,5, 9| , 
however  their  validation  with  experimental  data  has  not  often  been  attempted,  mainly  due  to  the  difficulty 
in  getting  experimental  results. 

The  present  work  includes  a theoretical  analysis  of  the  evolution  of  a water  spray,  given  an  initial 
measured  droplet  distribution,  characterised  by  a modified  Rink  distribution  function  |8|,  and  a comparison 
of  the  predicted  values  with  the  measured  ones  at  three  different  downstream  locations  along  the  spray  axis. 
The  droplet  distribution  has  been  divided  into  a discrete  number  of  size  groups  with  single  droplet 
equations  applied  to  each  one  of  them.  The  droplet  size  measurements  were  made  using  an  optical  technique, 
based  on  the  diffraction  pattern  formed  when  the  droplets  are  illuminated  by  a parallel  beam  of 
monochromatic,  coherent  light  | 7 ,9 , 11 | . 


2.  THEORY 

The  assumptions  used  where  the  following 
1:  Droplet  with  spherical  symmetry 


2:  Large  distance  between  droplets 

3:  Physical  properties  varying  with  temperature 

4:  Constant  pressure 

5:  Uniform  droplet  temperature 

6:  No  radiation  heat  transfer 


2.1  Droplet  heat  balance 

The  heat  balance  of  a droplet  exposed  to  a hot  environment  is  expressed  by  the  equation: 


h A (T  -T.)  - M.C  . 
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The  convected  heat  is  equal  to  the  increase  of  internal  energy  plus  the  latent  heat  of  vaporization. 


In  the  first  period  of  the  droplet  life  the  vaporization  rate  dM/dt  is  small  compared  with  the  increase 
in  internal  energy,  so  it  can  be  assumed  that  the  convected  heat  is  equal  to  the  increase  in  internal 
energy  of  the  droplet.  When  the  droplet  temperature  reaches  the  equilibrium  vaporization  temperature,  the 
temperature  remains  constant  (dTd/dt-0)  and  the  rate  of  heat  transmission  determine  the  vaporization  rate 
of  the  droplet.  As  the  air  temperature  is  much  greater  than  the  droplet  boilling  temperature,  the 
equilibrium  temperature  is  the  boilling  temperature  at  ambient  pressure. 


From  equation  (1),  the  droplet  heating  time  is  given  by: 
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Tdo  - Droplet  temperature  at  the  beginning 
Td  - Droplet  temperature  at  the  instant  t 

and  the  evaporation  rate  is 


2n(l + B) 


(2) 


(3) 


C (T  -T  ) 

Where  B is  Spalding  mass  transfer  number  defined  as  p^  °°  d 

h£g 

For  the  evaluation  of  the  heat  transmission  coefficient  h the  empirical  expression  from 
Marshall,  which  applies  a convection  correction,  was  used 


Ranz  and 


Nu«“  - 2 + 0,6  Re1/Z  Pr2/3 


(4) 


The  heat  transmission,  to  the  droplet,  is  a function  of  the  Reynolds  number,  it  increases  with  the 
velocity  as  well  as  with  the  droplet  diameter,  so  a drag  factor  must  be  evaluated  and  used  to  calculate 
the  change  of  the  droplet  Reynolds  number  with  time. 


2.2  Droplet  equation  of  motion 

The  droplet  motion  equation  is: 


d(Mu)  u2 

dt  "P8  2 


TT  R 


C.  + R 
d s 


(5) 


Rg-  source  term,  includes  the  gravity  forces  and  others , neglected  in  this  case 
Cj-  is  the  drag  coefficient,  evaluated  from  the  empirical  relation  proposed  by  Mellor 

Cd-  28/Re0,85  + 0.48. 

To  determine  the  relative  velocity  between  the  droplet  and  the  gas  the  initial  conditions  must  be 
known.  The  air  velocity  can  be  evaluated  because  the  mass  flow  and  the  temperature  are  measured.  The 
droplet  velocity  at  the  beginning  is  not  measured,  but  the  atomizer  used  in  the  experiments  was  a twin- 
-fluid  atomizer,  described  in  |8|,  and  the  initial  droplet  velocity  was  assumed  to  be  153!  of  the  mean 
stream  velocity.  As  equation  (5)  shows,  the  droplet  velocity  in  the  control  volume  is  a function  of  the 
droplet  diameter,  so  in  the  measuring  volume  the  droplet  concentration  is  a function  of  the  distance  from 
the  injector.  Figure  1 shows  the  droplet  volume  distributions  in  a spray  evaluated  at  four  different  axial 
stations.  In  figure  2 the  evolution  of  the  relative  droplet  velocity,  for  different  droplet  diameters,  can 
be  seen. 


2.3  Diffracted  Pattern 


With  the  equations  of  motion  and  heat  balance  it  is  possible  to  evaluate  the  droplet  size  distribution 
in  the  control  volume.  Knowing  this  distribution  the  evaluation  of  the  diffracted  light  pattern  from  a 
coherent  light  beam  crossing  the  control  volume,  is  made  using  the  Fraunhofer  theory  |ll|.  The  light 
diffracted  by  the  smallest  droplets  (D<20  pm),  for  which  this  theory  can  not  be  used,  is  evaluated  using 
the  Mie-theory. 


16-3 


4 

! 


3.  EXPERIMENTAL  SET  UP  AND  MEASUREMENTS 

Figure  3 shows  the  experimental  set  up.  In  a rectangular  channel  (20x  30  mm)  an  air  stream  flows  at 
controlled  temperature  and  velocity.  In  this  flow  a water  film  is  injected  and  atomization  occurs  by  the 
shearing  action  of  the  high  velocity  air  stream  against  the  low  velocity  liquid  film.  The  length  of  the 
channel  can  be  changed  to  allow  the  measurement  of  the  droplet  distribution  and  concentration  at  four 
different  downstream  distances,  10,  30.  70,  110  mm.  The  droplet  distribution  and  concentration  is  measured 
with  an  optical  method  described  in  j 7 1 . 

A parallel  monochromatic  light  beam  10  mm  diameter  from  a 5 mW  He.Ne  laser  system,  crosses  the  spray. 
Some  of  the  light  is  diffracted  due  to  the  droplets,  the  diffraction  pattern  depending  on  the  droplet 
size.  A converging  lens  focuses  the  non  scattered  ligth  in  the  centre  of  a receiving  plane,  where  a 
photomultiplier  measures  the  light  intensity.  Measuring  this  value  with  and  without  spray,  enables  the 
evaluation  |11|  of  the  liquid  volume  concentration  Cv,  based  on  the  transmission  law  for  polydispersed 
particles. 


E/E0  = exp(  - 3/2K  Cv  Jl/D32> 

K - Mean  scattering  coefficient 

l - Spray  dimension,  crossed  by  the  laser  beam  (see  fig.  3) 

Around  the  focus  of  the  same  converging  lens  a pattern  of  rings  of  the  diffracted  ligth  is  formed. 

Nine  concentric  circular  sectors  of  photomultipliers  measure  the  diffracted  ligth  intensity. 

Ten  sample-and-hold  units  were  used  to  garantee  that  an  instantaneous  light  profile  from  the 
photomultipliers  was  obtained.  The  signals  were  digitized  with  an  A/D  converter  and  the  data  was  stored 
in  paper  tape  for  later  processing  in  a BS  2000  (Siemens)  digital  computer. 

With  this  system,  the  diffracted  light,  which  has  been  shown  to  be  a function  of  the  droplet  size 
distribution,  and  the  transmitted  light,  which  is  related  to  the  liquid  concentration  in  the  spray,  were 
both  measured. 

3.2  Measurements 

Three  sets  of  measurements  have  been  made,  the  first  one  at  ambient  temperature,  one  at  200°C  and  one 
at  400°C,  all  at  atmospheric  pressure. 

The  influence  of  the  air  velocity  and  water  flow  rate  ir.  the  quality  of  the  spray  was  first  studied. 

Figure  4 shows  that  increasing  the  air  velocity,  the  spray  Sauter  Mean  Diameter  decreases,  and 
increasing  the  water  flow  rate  the  Sauter  Mean  Diameter  increases,  but  only  for  the  smaller  air 
velocities  this  parameter  is  important.  With  air  velocity  of  120  m/s  the  change  in  the  water  flow  fr<  a 5 
to  20  1/h  has  no  influence  in  the  S.M.D. 

These  results  are  in  good  agreement  with  the  results  of  Mullinger  and  Chigier  (10). 

The  air  velocity  is  also  important  in  the  liquid  volume  concentration  as  figure  5 shows.  Increasing 
the  air  velocity  the  liquid  volume  concentration,  measured  at  the  same  distance  (30  mm)  from  the  injector, 
decreases.  The  liquid  concentration  is  proportional  to  the  water  flow  rate  and  inversely  proportional  to 
the  air  flow. 

Figure  6 shows  the  influence  of  the  distance  in  the  Sauter  Mean  Diameter  and  in  the  liquid  volume 
concentration  for  an  air  velocity  of  120  m/s  and  ambient  temperature.  The  S.M.D.  remains  constant  and  the 
liquid  volume  concentration  decreases.  The  same  can  be  observed  in  figures  7 and  8 for  two  other  cases, 
air  velocity  180  and  190  m/s  and  air  temperature  400  and  200  C respectively  and  water  flow  rates  of  10, 
15,  20  1/h.  In  figure  8 the  predicted  change  in  liquid  concentration  with  distance  from  the  injector  is 
also  shown  (solid  lines).  The  predicted  evolution  is  based  on  the  first  measurement  of  each  curve. 

The  liquid  concentration  in  the  control  volume  changes  with  axial  distance  and  with  the  air  flow 
parameters.  This  is  due  to  the  different  droplet  velocities,  the  drag  force  changes  with  the  air  velocity 
and  temperature  and  with  the  droplet  diameter. 

As  figure  7 shows  the  droplet  diameter  does  not  change  with  axial  distance,  no  evaporation  or 
coalescense  occurs.  Predictions  of  the  variation  of  the  droplet  diameter  with  distance  from  injector, 
are  shown  in  figure  9,  evaluated  for  case  B. 

We  can  see  that  although  case  B and  C have  almost  the  same  velocity  the  S.M.D.  in  case  C is  smaller, 
this  being  due  to  the  higher  air  temperature  in  case  B that  causes  a worse  atomization.  For  the  same 
reason  cases  A and  B have  the  same  S.M.D.  although  the  air  velocities  are  quite  different. 

Figure  10  shows  the  measured  and  the  predicted  normalised  light  profiles  for  the  three  different 
experiments  (A,  B,  C),  at  30  mm  from  the  injector. 


Table  1 

Experimental  conditions 


Case 

A 

B 

C 

Air  vel.  (ra/s) 

120 

180 

190 

Air  temp.  (°C) 

20 

400 

200 

d32 

45 

43 

31 

Table  1 defines  the  parameters  of  experiments  A,  B,  C and  the  measured  Sauter  Mean  Diameter. 
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Table  2 contains  the  measured  and  predicted  values  of  liquid  volume  in  the  control  volume. 

Table  2 

Liquid  volume  variation  in  the  control  volume  (Z) 

Water  flow  rate  10  2/h 


dist.  (mm) 

A 

B 

C 

10 

100 

100 

100 

30 

81 

66(74) 

72(72) 

70 

67 

50(47) 

53(54) 

110 

61 

42(44) 

44(42) 

The  experimental  results,  (in  brackets)  for  a water  flow  of  10  1/h  are  in  good  agreement  with  the 
theoretical  ones. 

Case  B is  where  the  deviation  from  predictions  is  greater  (see  fig.  8)  due  to  the  higher  temperature  of 
the  air  stream  that  causes  a shear  layer  of  two  different  mediums  to  interfere  with  the  light  beam. 

In  some  cases,  interference  of  the  oscillating  shear  layer  between  the  hot  and  cold  air  was  noticed. 

3.3  Source  of  errors 


As  figure  10  shows,  in  the  first  inner  ring  the  light  intensity  does  not  agree  with  the  predictions. 
This  is  due  to  the  small  radius  of  the  inner  ring  and  the  presence  of  some  scattered  light  even  during 
calibration,  made  without  spray. 

We  have  checked  the  influence  of  the  oscilating  interface  between  the  hot  and  cold  air,  and  it  has 
been  noticed  that  this  influence  principaly  disturbs  the  first  inner  ring.  It  has  also  been  noticed  that 
this  influence  can  be  minimised  with  glass  windows  in  the  air  channel. 

In  all  the  measurements,  the  first  inner  ringhave  not  been  taken  into  account. 


4.  CONCLUSIONS 

Measurements  of  the  Sauter  Mean  Diameter  of  a water  spray  from  a twin-fluid  atomizer  have  been  made, 
for  different  air  temperatures  and  velocities,  using  an  optical  method. 

The  influence  of  the  air  velocity  and  temperature  and  of  the  water  flow  rate  in  the  spray  S.M.D.  and 
in  the  liquid  volume  concentration  at  different  distances  from  the  injector  have  been  studied. 

The  results  agree  with  the  predicted  ones  and  with  previous  measurements  made  in  similar  atomizers 

(10,  12). 

Care  must  be  taken  when  the  light  beam  passes  across  regions  of  varying  density,  as  in  the  case  of 
heated  air  streams  or  flames. 
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Figure  i.  Droplet  volume  distribution  evaluated 
at  four  axial  stations. 
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DISCUSSION 


P.Ramette,  Fr 

Vous  avez  prdsente  des  correlations  en  fonction  de  la  vitesse  de  1’air  et  du  debit  du  jet.  Avez-vous  essaye  d’exprimer 
vos  rtsultats  en  fonction  de  la  vitesse  relative  du  jet  par  rapport  a Pair? 

Rdponse  d’Auteur 

Non.  Nous  avons  choisi  comme  variables  cedes  qui  correspondent  aux  mesures  directes. 


